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, Abstract. 

' We study backreaction analytically using the parabolic Lemaitre-Tolman-Bondi 

1} ■ universe as a toy model. We calculate the average expansion rate and energy density 

Q \ on two different hypersurfaces and compare the results. We also consider the Hubble 

law and find that backreaction slows down the expansion if measured with proper time, 
but speeds it up if measured with energy density. 
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The homogeneous and isotropic Friedmann-Robertson- Walker models are usually 
thought to give a good description of the average behaviour of the universe, since the 



1. Introduction 



tjr)! universe appears to be homogeneous and isotropic on sufficiently large scales (though see 
PP). The reasoning behind the FRW equations is that one takes the average of the real 
(inhomogeneous and anisotropic) metric and energy-momentum tensor, and plugs these 
into the Einstein equation. However, the physically correct thing would be to first plug 
the real inhomogeneous and anisotropic quantities into the Einstein equation, and only 
then take the average. Because the Einstein equation is non-linear, these two procedures 
are not equivalent. In other words, the averages of the real quantities do not satisfy 
the Einstein equation. The feature that the average behaviour of an inhomogeneous 
spacetime is not the same as the behaviour of the corresponding smooth spacetime 
(that is, one with the same average initial conditions) is called backreaction. In the 
context of cosmology, the issue was highlighted in [2j (it had been discussed earlier in 
[3]) as the fitting problem: how do we find the homogeneous and isotropic model which 
best fits the real inhomogeneous and anisotropic universe? 

The relation between the average sources (energy density, pressure, etc.) and 
average geometric quantities (expansion rate, shear, etc.) in a general spacetime filled 
with a perfect irrotational fluid is known |1J E] . However, the system of equations is 
not closed, which means that different inhomogeneous spacetimes with the same initial 
averages evolve differently even as far as average quantities are concerned. Since there 
is no procedure for finding the average behaviour of a given spacetime (short of solving 
it exactly), it is difficult to quantitatively evaluate the importance of backreaction in 
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cosmology. One notable calculation is jH], where the behaviour of an universe-sized box 
filled with inhomogeneous expanding dust is followed numerically using the Newtonian 
limit of the exact backreaction formalism of Hj. The results show unambiguously 
that backreaction is a real phenomenon which can have a large impact: for example, 
it can make regions with no initial overdensity turn around and collapse. However, 
the Newtonian limit is not expected to fully capture the relativistic backreaction. (In 
particular, the global backreaction vanishes identically for periodic boundary conditions, 
unlike in the relativistic case.) 

Fully relativistic quantitative studies of backreaction have usually been done in 
the context of perturbative solutions around a FRW background jSl El EH UH EE21 EH 
El EH EE E] (see J3| for a more extensive list of references). Of particular interest 
has been the possibility that backreaction of long wavelength perturbations could lead 
to a dynamical relaxation of vacuum energy, at the same stroke providing an elegant 
inflationary mechanism and explaining why the vacuum energy is so much smaller than 
theoretically expected [T31 [T3 [THJ UH1 EOJ ■ It has also been suggested that this could solve 
the coincidence problem of the effective vacuum energy density being of the order of the 
matter energy density today [THI ED] • Another possible explanation for the coincidence 
problem is backreaction from structure formation, that is, from small wavelength modes 
[HI EE EH 121] • Since the acceleration of the universe seems to have started around the era 
when structure formation is important, it seems a natural possibility that this deviation 
from the simple prediction of deceleration in the homogeneous and isotropic models with 
normal matter could be related to the growth of inhomogeneities in the universe. 

Perturbative treatment of backreaction has two drawbacks. First, naive 
perturbative results break down when perturbations have an effect on the background, 
in other words when backreaction is important. In addition to the usual issues of 
cosmological perturbation theory such as gauge-invariance |H1 El 122] j one has to worry 
about new problems such as convergence and consistency of the perturbative expansion. 
These make consistent backreaction calculations in the perturbative framework an 
involved task. Second, present-day universe is not perturbatively close to homogeneity, 
but contains non-linear structures. If the coincidence problem is to be solved by 
backreaction from structure formation, one would have to go beyond perturbation theory 
(though if the effect involves mainly perturbations breaking away from the linear regime, 
a quasi-perturbative treatment might still be possible). 

We will study backreaction analytically and free of perturbative ambiguities with 
an exact toy model, the Lemaitre-Tolman-Bondi model [23 123 I2S] (see also [211 1211) ■ 
The LTB model is the spherically symmetric dust solution of the Einstein equation 
(or rather the family of such solutions). Like the perturbed FRW spacetime, the LTB 
spacetime is a generalisation of the FRW universe. The LTB model can be viewed as 
an Einstein-de Sitter universe with a single spherically symmetric perturbation which 
can be arbitrarily large, as opposed to the linearly perturbed FRW universe which has 
an ensemble of small perturbations. The LTB model can describe the collapse of an 
overdensity or the formation of a void in an expanding universe |2*5j . and it has also 
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been used to describe the local inhomogeneous universe [29] . As a model for the entire 
universe, it is a toy model which is useful because the backreaction problem can be 
studied quantitatively and without any approximations, as has been previously done in 
[3Ti] (averaging in the LTB model has also been discussed in [21]; see also [22])- 

In a perturbative framework one aspect of backreaction is that inhomogeneities 
change the behaviour of a smooth background. This feature is obviously not present in 
an exact inhomogeneous solution, but the difference between the average behaviour and 
the behaviour of the corresponding smooth spacetime can be studied unambiguously. 
With an exact solution one can also study the choice of the hypersurface of averaging 
in isolation of the different issue of the choice of gauge [HI E3 UH U21 EH 122] • 

The present study is a simple exact non-perturbative counterpart of [12] , where we 
calculated perturbatively the average expansion rate of a linearly perturbed Einstein-de 
Sitter universe. In section 2 we calculate the average expansion rate in the simplest 
LTB solution for two different choices of hypersurface, and compare the results with 
each other and the FRW universe. In section 3 we look at the impact of backreaction 
on the Hubble law and discuss our results. 



2. The backreaction calculation 

2.1. The Lemaitre-Tolman-Bondi model 

The metric. We will proceed as in ^2]: we will write down the metric, find the proper 
time and the expansion rate and take the average over the hypersurface of proper time. 
We will also take the average over the hypersurface of constant coordinate time, and 
compare the results. 

The LTB model is the most general spherically symmetric dust solution. The 
Einstein equation is (we take the cosmological constant to be zero) 

1 1 

where Mpi = 1/a/87tG/v is the (reduced) Planck mass, p is the energy density and 
is the velocity of the matter fluid, with = —1. The spherical symmetry allows the 
Einstein equation (Q) to be solved exactly. The metric turns out to be 

ds 2 = -dt 2 + ^it- dr 2 + R{t, r) 2 {d6 2 + sin 2 6dtf) , (2) 

where the functions R(t, r) and E(r) are related to each other and to the energy density 
p(t, r) as follows 

no 1 m(r) , , 

/, \ rn'jr) 
P(t ' r) =R(t,ry R %r) ' (3) 
where dots and primes denote derivatives with respect to t and r, respectively, and 
m(r) is a function which describes how much energy there is within the radius r. (The 
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freedom to choose the r-coordinate means that m(r) can be redefined at will.) The 
velocity of the matter fluid is simply 



There are three different classes of solutions to Q, for E > 0, E < and E = 0. 
We will consider the solution with E = 0, known as the parabolic solution, which is the 
simplest and the most analogous to the Einstein-de Sitter universe. We assume that 
ml > and choose the radial coordinate such that m(r) = AMp t r 3 j (9tf) , where t\ is a 
positive constant with the dimension of time. We then have 



where t (r) is an arbitrary function. It is transparent that the solution reduces to the 
FRW case if to(r) = constant. There is a singularity at t — to(r), and we have chosen 
to look at the case t > t (r), which describes matter expanding away from the big bang 
rather than collapsing towards a future singularity. The function to(r) is the big bang 
time for a comoving observer at constant r: different observers find their part of the 
universe to have emerged from the singularity at different times in the past. In order to 
avoid a singularity caused by shells of matter crossing, we must have t' (r) < [33] . 

The proper time. As emphasised in JTQlEIlGII], it is important to cast things in terms 
of the proper time of the observer. Given (jlj), the derivative in the direction orthogonal 
to the hypersurface defined by the velocity of comoving observers is d T = m m V m = W 
From the condition d T r = 1 we get the proper time r = t + f(r), where f(r) is an 
arbitrary function. 

In the perturbed FRW case [13] there was a difference between the background 
coordinate time t and the proper time r because the velocity of the perturbed 
comoving observers was different from the background velocity, leading to a perturbed 
t. Here the situation is different. The velocity of comoving observers is the same as 
in the smooth case, but different observers start their clocks at different times. When 
discussing the expansion rate measured by cosmological observers, we want to compare 
observers whose clocks show the same time. We therefore choose f(r) = —t (r), so the 
proper time is 



which is the "local age of the universe", the time since the big bang measured by a 
comoving observer at constant r. Note that in order for r to be a time coordinate, the 
normal V M r has to be timelike, which translates into the constraint R'(t,r) > —t' (r). 

The expansion rate. We are interested in the expansion rate measured by a comoving 
observer, given by 



(1,0,0,0) . 



(4) 




(5) 



T (t,r) = t-t (r) , 



(6) 



0(*,r) = VX» 



(7) 
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which in the FRW case reduces to 9 = 3H, where H is the Hubble parameter. Plugging 
in (0), and © we have 

^r) = Ul--^—) . (8) 



t \ 3r + 2rr' / 

The expansion rate of an Einstein-de Sitter universe is 2/r, so the second term is 
the contribution of the inhomogeneity. Since r' = —t' > 0, this term is negative, and 
the expansion in terms of the proper time is always slower than in the FRW case, and 
slows down with increasing r. In the converse situation with r' < 0, the expansion rate 
would be larger than in the FRW case, but the inner shells of matter would overtake 
the outer shells, resulting in a singularity. (In a more realistic model, shell-crossing 
would mean that the decription of matter as a pressureless ideal fluid breaks down, not 
necessarily that there is a singularity.) 

The energy density. From (J3J) and (jSJ) we have 
, , AM 2 3r 

In the solution with E = 0, the density is uniquely determined by the local big 
bang time to, and one can see from p concretely what kind of a physical situation a 
given to corresponds to. In particular, when t is constant, Q reduces to the FRW 
case, p = 4Mp;/(3r 2 ). It is possible to specify a LTB solution by giving the energy 
density on initial and final hypersurfaces [2E], instead of giving E and t as done here. 



2.2. Taking the average 

The average expansion rate. To obtain the average expansion rate, we should integrate 
the local expansion rate (jHJ) over the appropriate hypersurface. To find the integration 
measure on the hypersurface of constant r, let us write the metric (J2J) in terms of r: 

ds 2 = - dr 2 - 2t drdr + (R> 2 - t' 2 )dr 2 + R 2 (d6 2 + sin 2 a <*^ 



= - ^fr^^ 2 + ( R ' 2 ~ f o) ( dr - w^ dr ) 2 + R2 ^° 2 + sin2 9d ^ ' ( 10 ) 

where we again see that R' > —t' is required for r to be a time coordinate. We can now 
read off the integration measure on the hypersurface of constant r [SI] : 

|det^| 1/2 = (R' 2 -t' 2 ) l ^R 2 sin 9 . (11) 
The average of a scalar observable F(t,r) over the hypersurface of constant r is 

then 

(F) Tc 



fd 4 2 


det ^g i3 


^6(t-t {r)-r c )F{t,r) 


Jd 4 x|det ( T 


9ij 


V 2 5(t - f (r) - r c ) 



J °° drr 2 



Toy I 3T c -2rf Q \ _ t , 2 



1/2 

F(t (r) + r c ,r) 

— , (12) 



1/2 
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where 8(t — to(r) ~ T c) is the delta function, r c is the constant value of proper time which 
labels the hypersurface, and on the second line we have plugged in © and (fTTjl. To keep 
the analogy with the FRW case as close as possible, we will consider only spacetimes 
where the coordinate r ranges from to oo. 

In comparison, the integration measure on the hypersurface of constant t is 

|det ( %-| 1/2 = R'B 2 sin 6 , (13) 

and the average over the hypersurface of constant t reads 

f d 4 x| det %/ /2 5(t - t e )F(t,r) 



(F)t 



f d 4 :r|det ^g l3 \ 1 ' 2 5(t - t c ) 
r °° drr 2 (t c - t )[3(t c - t ) - 2rt' )F{t c , r) 



j °° drr2(t c - t )[3(t c - t ) - 2rt ] ' (U) 

where t c is again the constant labelling the hypersurface. There are two differences 
between the averages (|T2*j) and (|14|). The integration measure is different, and the 
quantity being held constant in the integral is different. In general, both are expected 
to be important, but in the present case the integration measure will in fact turn out 
not to matter. We will now evaluate the average of the expansion rate 6 over both 
hypersurfaces for some simple choices of t . 



Comparing the averages. Our LTB solution is specified once we give the function to, 
which has to satisfy R' > —t' > 0. We consider three examples of to- They are given in 
Table ^ along with the average expansion rates and average densities (tq is a positive 
constant). The exponential factors in the first two cases are needed for satisfying the 
condition R' > — 1 ; without them, the averages would be the same but the hypersurfaces 
would not be spacelike for all r c . We also give the effective scale factor a T defined by 
3d T a T /a T = (9) T , in analogy with |H El EE E|i (h is a positive constant). 



Table 1. Average expansion rates and densities. 





to(r) 


(0)r 




(P)r 


(P)t 


a T 




Case 1 


—re ~ 


2 1 
T 2 











(r/h, 


jl/3 


Case 2 
Case 3 


r e-^ln r+r ° 

ro 


2 3r+r 
t 3r+2ro 

2 




2 


3t 2 3r+2r 




AMI, 


[r(r- 

(r/h] 


h2r /3)/tI]V3 

)2/3 


l+r 3 /{2t(ro) 


T 


f 


3r 2 


3t 2 



In cases 1 and 2 the big bang occurred at t = in the center at r = 0, while the 
far away shells of constant r asymptotically approach being infinitely old. In contrast, 
in case 3 there is only a finite difference (of ro) in the big bang time between the center 

| Note that this definition is different from the one used in [T3HT?)] . 
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and the region at asymptotic infinity. The three cases cover the possible asymptotic 
behaviours for the quantity rr'\ divergent, finite and zero. 

Scalar quantities (such as the expansion rate) are invariant under coordinate 
transformations, but spatial averaging is not a covariant procedure, and the choice of 
hypersurface makes a physical difference. Cases 1 and 2 illustrate the drastic difference 
between the hypersurfaces of constant r and constant t. In case 1 the average expansion 
rate on the hypersurface of constant r is 1/2 of the FRW rate, and in case 2 it starts 
at 1/2 of the FRW rate and asymptotically approaches it as r grows, and the scale 
factor correspondingly interpolates between (r/ti) 1 ' 3 (case 1 ) and (r/tx) 2 / 3 (the FRW 
case). However, in both cases 1 and 2 the average expansion rate on the hypersurface of 
constant t is zero. In case 3 the average over both hypersurfaces yields the unperturbed 
expansion rate in terms of the appropriate time. 

It may seem strange that the average of the expansion rate is zero, given that the 
expansion rate is positive definite. The resolution is that the average is dominated by the 
region at asymptotic infinity, since it contains infinitely more observers than the region 
within any given finite radius. Since the age of the region at large r is asymptotically 
infinite, its expansion rate is asymptotically zero, so the average expansion rate goes to 
zero. It is for the same reason that the averages turn out to be the unperturbed values 
in case 3: the big bang time to goes asymptotically to zero, so there is asymptotically no 
difference between r and t, and therefore no backreaction. In fact, the averages over the 
two hypersurfaces are simply the asymptotic limits of the quantities (JBJ), © with the 
appropriate time coordinate held fixed. Since the three cases we have studied cover all 
possible limits of the quantity rr', we have exhausted the possibilities for the averages: 
any permissible function to leads to one of the average quantities in cases 1 to 3 (apart 
from a trivial global change in the normalisation of t). The averages over all space do 
not depend on the integration measure, since in the asymptotic limit the measure in 
the nominator and the denominator cancels. This feature is peculiar to the spherical 
symmetry of the LTB model. In a realistic model of the universe where all points are 
statistically equivalent, the integration measure is expected to make a difference. It 
would also be more physical to average over the finite region that has been in causal 
contact with a given observer, in which case the backreaction would not depend only 
on the asymptotic limit, and would be finite in all cases. 

The above examples show that even in the LTB model, where (unlike in the 
perturbed FRW universe) the coordinate time t is a physically measurable quantity, 
taking averages over the hypersurface of constant t can give misleading results. We are 
interested in the average expansion rate measured by local observers as a function of the 
proper time they measure, so we should compare observers whose clocks show the same 
time. Assuming that cosmological observers normalise their clocks to the locally inferred 
big bang time to instead of synchronising them to some global constant, observers further 
out in r will reach the same value of r earlier in terms of t: the averaging surface tilts 
to the past with increasing r. 
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3. Discussion 

The Hubble law. To get another viewpoint on the backreaction, let us look at the 
relation between the average expansion rate and average density in case 2. From the 
expressions given in Tabled we can solve for r in terms of {p) T , and insert this into the 
expression for {9) T . The resulting Hubble law is 



where we have written the square of the expansion rate in three different forms: as a 
function of the proper time, the average energy density and the effective scale factor. 
The FRW relations are recovered as vq/t goes to zero. 

From the first expression, (|15|) . it seems that the expansion is slower than in the 
homogeneous case. However, from the second and third expressions, (j!6|) and (fTTj) . it 
seems as if the universe would instead expand faster as there is an additional positive 
(p)^-term driving the expansion. We have the paradoxical-seeming situation that the 
expansion is slower than in the FRW case if measured in terms of the proper time, but 
faster if measured in terms of the energy density (or the scale factor, since conservation 
of mass implies (p) T oc a~ 3 )- The reason is that the backreaction reduces (in terms of r) 
(p) T more than it reduces (9) T , as we see from TableHJ Therefore, the dependence of (9) T 
on (p) T has to be stronger than in the FRW case to get the dependence on r right. This 
underlines the point, mentioned in [TIHll3j. that one should be careful to identify which 
quantity is being measured in terms of what parameters when considering backreaction. 

A measure of backreaction. In jT3j it was suggested that since the Weyl tensor gives 
a measure of the inhomogeneity and anisotropy of a spacetime it could provide some 
indication of whether backreaction is important. For the LTB solution with E = 0, the 
ratio of the square of the Weyl tensor to the square of the scalar curvature is, from (J2J) 
and ©, 



which does measure the importance of backreaction. For example, in case 2 it is, when 
averaged over the hypersurface of constant r, essentially the ratio (r /r) 2 . Note that 
the information entropy measure introduced in [33] to quantify the inhomogeneity of a 
spacetime measures the difference between a spacetime and its average, not the difference 
between the average and the corresponding smooth spacetime. (At any rate, the measure 
was defined for compact domains, and so is not applicable to the present case.) 

One can calculate the backreaction for geometric quantities other than the 
expansion rate and the Weyl tensor. One interesting candidate is the shear, which 




(15) 



(16) 



(17) 




(18) 
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is zero in the FRW case. Like the square of the Weyl tensor, the shear turns out to 
be a measure of the backreaction both in the LTB case and in the perturbative case of 
[T3] . In the perturbative case, the shear becomes of the order of the scalar curvature, 
just like the Weyl tensor. This is in conflict with observations, but the perturbative 
analysis of ^3] is at any rate expected to break down when backreaction is indicated to 
be important. 

Conclusion. We have calculated the average expansion rate and energy density in the 
simplest, parabolic, LTB solution analytically for some example cases (which turn out 
to cover all possibilities). Comparison of the averages with each other demonstrates the 
importance of the choice of hypersurface, and comparison with the FRW case shows 
that backreaction slows down the expansion if measured in terms of the proper time, 
but speeds it up if measured in terms of the energy density or the scale factor. The 
calculation is an example of an exact, quantitative study in backreaction. 

Unlike the perturbed FRW metric, the LTB metric is not meant to be a realistic 
model of the universe, but it is a useful toy model for studying backreaction because one 
can obtain analytical results. One way towards a more realistic backreaction calculation 
from the first order perturbed FRW model used in J3| is to do a consistent second order 
perturbative analysis. Another possibility could be a quasi-perturbative approach where 
one keeps to first order perturbation theory for scales which are in the linear regime, 
but applies a non-perturbative solution for scales which have gone non-linear. Since 
the LTB model can describe the formation of clusters and voids smoothly starting from 
small perturbations j2H], it is ideal for this purpose. (A simple version of embedding 
LTB solutions into a FRW universe was done in [HE], and the idea goes back to |23j-) 
This would involve the LTB solutions with E ^ 0, and it would be interesting to first 
study how the average expansion rate behaves in those solutions themselves. 
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